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Résumé Le problème du conflit intrinsèque à la combinaison d'informations a poussé à de 
nombreuses réflexions ces dernières années, en particulier dans le cadre de la théorie des fonc- 
tions de croyance. Nous pouvons résumer les solutions apportées par trois façons de considérer 
le problème : premièrement, nous pouvons chercher à réduire voire supprimer le conflit avant 
la combinaison d'informations, deuxièmement nous pouvons gérer le conflit de façon à ce qu'il 
n'intervienne pas lors de la combinaison et nous le considérons seulement lors de la décision 
finale, et troisièmement nous pouvons tenir compte du conflit lors de l'étape de combinaison. 
Si la première solution paraît la meilleure elle n'est pas toujours réalisable ou suffisante. Il peut 
être difficile de chercher à départager philosophiquement les deux dernières stratégies. Cepen- 
dant d'un point de vue applicatif seule la décision compte, et c'est donc dans cette optique que 
nous chercherons à les comparer. Nous proposons ici une nouvelle règle qui a pour principe de 
répartir le conflit proportionnellement sur les éléments produisant ce conflit. Nous comparons 
les différentes règles à partir de données réelles en imagerie Sonar et en classification de cibles 
Radar. 

Mots Clés : Fonctions de croyance, Conflit, Règle de Combinaison, Imagerie 
Sonar, Classification de cibles Radar 

> ■ 

^ ■ Abstract Thèse last years, there were many studies on the problem of the conflict coming 
from information combination, especially in évidence theory. We can summarise the solutions 
for manage the conflict into three différent approaches : first, we can try to suppress or reduce 
the conflict before the combination step, secondly, we can manage the conflict in order to give no 
influence of the conflict in the combination step, and then take into account the conflict in the 
décision step, thirdly, we can take into account the conflict in the combination step. The first 
approach is certainly the better, but not always feasible. It is difficult to say which approach 
is the best between the second and the third. However, the most important is the produced 
results in applications. We propose here a new combination rule that distributes the conflict 
proportionally on the élément given this conflict. We compare thèse différent combination rules 
on real data in Sonar imagery and Radar target classification. 

Key Words : Belief Functions, Conflict, Combination Rule, Sonar Imagery, Ra- 
dar Target Classification 
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1 Introduction 



La fusion d'informations a depuis plusieurs années permis d'apporter des solutions à la 
combinaison d'informations issues de diverses sources afin d'améliorer la prise de décision par 
exemple dans un système de classification. Le fait de devoir tenir compte de plusieurs sources a 
poussé le développement de plusieurs théories permettant de modéliser finement les informations 
issues de ces sources en terme de fiabilité, incertitude, imprécision ou encore incomplétude. 
Parmi ces théories de l'incertain, la théorie des fonctions de croyance issue des recherches de 
Dempster et Shafer pQCLZ] a connu ces dernières années un développement important notamment 
dans sa mise en forme pour de nombreuses applications. 

La théorie des fonctions de croyance est fondée sur la manipulation des fonctions de masse 
(ou masse élémentaire de croyance). Les fonctions de masse sont définies sur l'ensemble de 
toutes les disjonctions du cadre de discernement 6 = {Ci, . . . , C n } et à valeurs dans [0, 1], où 
les Ci représentent les hypothèses supposées exhaustives et exclusives. Cet ensemble est noté 
2 e . Généralement, il est ajouté une condition de normalité, donnée par : 

E = w 

où rrij(.) représente la fonction de masse pour une source (ou un expert) Sj, j = 1, ...,M. Les 
éléments X tels que m(X) > sont appelés les éléments focaux. La première difficulté est donc 
de définir ces fonctions de masse selon le problème. A partir de ces fonctions de masse, d'autres 
fonctions de croyance peuvent être définies, telles que les fonctions de crédibilité, représentant 
l'intensité que toutes les sources croient que les éléments focaux d'une source affirment la 
croyance en un élément. Elles sont données pour tout X G 2 e par : 

bel(X)= £ m(Y), (2) 

Y&2 X ,Y^V 

ou encore les fonctions de plausibilité, représentant l'intensité avec laquelle on ne doute pas en 
un élément, données pour tout X G 2 e par : 

pl(X) = E m(F) =bel(6)-bel(X c ) = l-m(0)-bel(X c ), (3) 

Y£2 e ,YnX^ê 

où X e est le complémentaire de X. 

Afin de conserver un maximum d'informations, il est préférable de rester à un niveau crédal 
(i.e. de manipuler des fonctions de croyance) pendant l'étape de combinaison des informations 
pour prendre la décision sur les fonctions de croyance issues de la combinaison. Si la décision 
prise par le maximum de crédibilité peut être trop pessimiste, la décision issue du maximum 
de plausibilité est bien souvent trop optimiste. Le maximum de la probabilité pignistique, 
introduite par [20], reste le compromis le plus employé. La probabilité pignistique est donnée 
pour tout X G 2 e , avec X ^ par : 



|xn Y| 


m{Y) 




y\ 


1 - m(0) 
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Notons que nous obtenons ainsi une probabilité peu conforme à la notion de fonction de masse. 
Le maximum de probabilité pignistique est généralement considéré uniquement sur les hypo- 
thèses singletons Ci à cause de l'additivité des probabilités. En effet : 



De la même façon la croissance des fonctions de crédibilité et de plausibilité pousse à prendre 
une décision uniquement sur les singletons. Si l'application autorise à prendre une décision finale 
sur une disjonction, il est alors préférable de considérer les fonctions de masse directement ou de 
développer une autre fonction de croyance non croissante adaptée à l'application. Dans un cas 
général une fonction de croyance non croissante sur le treillis d'inclusion mène à des résultats 
rapidement contre intuitif. Il peut également être intéressant de ne plus considérer l'unique 
hypothèse la plus vraisemblable, mais un sous-ensembles des hypothèses les plus vraisemblables 



Chercher à combiner des informations issues de plusieurs sources qui ne sont pas en accord 
fait bien souvent apparaître un conflit entre les sources. Ce conflit peut provenir d'un manque 
d'exhaustivité des sources, ou encore d'un manque de fiabilité de celles-ci [IJ. Plusieurs ap- 
proches sont envisageables pour gérer ce conflit. Dans le cas d'un manque d'exhaustivité des 
sources, une approche classique est la technique du hedging qui consiste à ajouter un élément au 
cadre de discernement. Lorsque la fiabilité des sources peut être estimée, une approche consiste 
à affaiblir les masses selon cette fiabilité. Cependant, depuis le problème posé par Zadeh, que 
nous reprendrons plus loin, de nombreuses règles de combinaison ont été proposées afin de 
manager le conflit [23l [71 [Hl ŒHl EU [121 ŒH El ŒEl E] - 

Nous proposons ici une nouvelle règle qui a pour principe de répartir le conflit propor- 
tionnellement sur les éléments produisant ce conflit. Dans un premier temps nous rappelons 
les différentes règles de combinaison, puis nous montrons en quoi il peut être intéressant de 
considérer cette nouvelle règle. Les différentes règles proposées au cours de ces dernières années 
peuvent être performantes selon les applications considérées. Dans une dernière partie, nous 
étudions le comportement de ces règles en terme de décision, dans le cadre de deux applications 
à partir de données réelles : en imagerie Sonar et en reconnaissance de cibles Radar. 

2 La combinaison 

Différentes approches de combinaison des fonctions de masse ont été proposées. La règle 
orthogonale de Dempster-Shafer non normalisée proposée par Smets pjj], est définie pour deux 
fonctions de masse vn\ et m 2 et pour tout A G 2 e par : 



betP(Ci U C 2 ) > betP(a), i = 1,2. 



M- 




(5) 



BnC=A 



Elle est donnée pour M experts, pour tout X £ 2 



par : 




lin...nr M =Xi=i 



(6) 
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où Yj G 2 e est un élément focal de la fonction de masse de l'expert j, et rrij(Yj) la masse 
associée. 

Cette opérateur est associatif et commutatif. La masse affectée sur l'ensemble vide s'inter- 
prète comme une mesure de conflit. Ce conflit peut provenir d'un manque d'exhaustivité des 
sources, ou encore d'un manque de fiabilité £[]. Dans le premier cas le fait d'avoir une masse 
non nulle sur l'ensemble vide (cas d'un monde ouvert) est concevable. Si l'on souhaite rester en 
monde fermé, la technique du hedging consiste à répartir le conflit sur le nouvel élément e, et 
la règle est alors donnée par : 

m h {X) = m Conj (X), VX^0 
m h {e) = m Co nj(0)- 

Si les sources ne sont pas fiables, lorsqu'il est possible de quantifier la fiabilité de chacune 
des sources il est important de procéder à un affaiblissement en redéfinissant les fonctions de 
masse par : 



m'jiX) = a jmj (X), VX G 2 e 

wf j (e) = i- aj (i- mj (e)). 



> ' J - ' (8) 



aj G [0, 1] est le coefficient d'affaiblissement de la source Sj qui est alors une estimation de la 
fiabilité de l'expert Sj, éventuellement comme une fonction de X G 2 e . 

Les experts pouvant s'exprimer dans 2 e (s'ils ne sont ni sûrs ni précis), l'apparition de 
conflit est inévitable. En effet, la règle orthogonale n'est pas idempotente, ainsi nous pouvons 
définir un auto-conflit d'ordre n pour chaque source j par : 

n 

a n {j) = © rnj(0), (9) 

k=l 

où l'opérateur © est l'opérateur de combinaison conjonctive de l'équation ((El) et nous avons la 
propriété : 

a n (j) < On+lO')- ( 10 ) 

Initialement, Dempster et Shafer ont proposé une règle orthogonale normalisée, afin de 
rester en monde fermé. Ainsi la répartition du conflit se fait de manière uniforme lors de la 
combinaison, pour tout X G 2 e , X ^ par : 

1 - ™Conj (0) y in ^ M=x f = \ 1 - ^Conj (0) 

où Yj G 2 e est la réponse de l'expert j, et rrijiYj) la fonction de masse associée. 

Smets ne répartit le conflit établi sur l'ensemble vide lors de la combinaison, qu'à l'étape 
de décision en prenant le maximum de la probabilité pignistique. Il multiplie toutes les masses 

par — r- (cf. équation ©). Ce critère de décision présente un compromis entre une 

1 - "'Co:,j (0) 

décision pessimiste par le maximum de crédibilité et une décision optimiste par le maximum de 
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plausibilité, aussi bien en monde ouvert qu'en monde fermé. Mais ces trois critères produisent 
la même décision que l'on normalise lors de l'étape de combinaison ou bien lors de la décision. 

La gestion du conflit peut cependant être réalisée lors de la combinaison de manière différente 
que dans le cas de la règle orthogonale normalisée, comme nous le verrons en rappelant l'exemple 
de Zadeh. Ainsi Yager [23J répartit le conflit sur l'ignorance totale (Le. sur la masse de 0) afin 
de rester en monde fermé et considère qu'on ne sait rien en cas de conflit. Ainsi la règle qu'il 
propose est donnée par : 

m Y (X) = m C onj(X),VAG2 e , X ^ 0, X^O 

m Y (e) = m Co nj(0) + m Conj (0) (12) 
m Y (0) = 0. 

Dans la plupart des applications, nous cherchons à prendre une décision sur les singletons 
et non sur l'ensemble 2 e . C'est ce qui est réalisé par la décision du maximum de probabilité 
pignistique. Ainsi, le problème soulevé par Zadeh (cf. tableau CD) lors d'un fort conflit n'est pas 
résolu, et la décision sur les singletons reste G. 



Tab. 1 - Exemple de Zadeh 








A 


B 


AU B 


c 





777.1 





0.9 








0.1 













0.9 





0.1 



















1 







0.99 











0.01 





my 














0.01 


0.99 



Dubois et Prade [7] ont proposé une gestion plus fine du conflit en répartissant le conflit 
partiel (par exemple issu uniquement de deux sources, l'une annonçant A et l'autre B) sur les 
ignorances partielles (c'est-à-dire AU B). Cette règle est donnée pour tout X G 2 e , X ^ par : 

M M 

m DP (X)= W m ^)+ E II" 1 ^). ( 13 ) 

Y x n...r\Y M =x j=i y 1 u...uy m = x i=i 

Y 1 n...nY M = 

où Yj G 2 e est un élément focal de l'expert j, et rrij{Yj) la fonction de masse associée. Si nous 
reprenons l'exemple de Zadeh, nous obtenons : 
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A 


B 


AU B 


G 


AuC 


BUC 





m DP 











0.99 


0.01 











bel DP 











0.99 


0.01 


0.01 


0.01 


1 


PW 





0.99 


0.99 


0.99 


0.01 


1 


1 


1 


betPDp 





0.495 


0.495 


0.99 


0.01 


0.505 


0.505 


1 



Ainsi avec le maximum de probabilité pignistique ou de la plausibilité nous choisissons A ou B 
sans pouvoir les distinguer à cause de l'égalité des masses données par les différents experts sur 
A et B. Si ces masses sont différentes la décision est alors possible, comme l'illustre l'exemple 
choisi suivant : 








A 


B 


AU B 


G 


AU C 


BUC 





777.1 





0.5421 


0.0924 





0.2953 








0.0702 


m 2 





0.2022 


0.0084 





0.6891 








0.1003 


m DP 





0.1782 


0.0106 


0.0233 


0.2815 


0.4333 


0.0662 


0.007 


bel D p 





0.1782 


0.0106 


0.2120 


0.2815 


0.8929 


0.3583 


1 


pW 





0.6417 


0.1071 


0.7185 


0.7880 


0.9894 


0.8218 


1 


betPop 





0.4088 


0.0577 


0.4665 


0.5335 


0.9423 


0.5912 


1 



Sur cet exemple la décision prise sera donc la classe C. 

De manière générale, la répartition du conflit peut s'écrire [TU! Œ21 ŒH] : 

m c (X) = mconj(X) + w(X)m Co nj(0), (14) 
où w(X) = 1. Toutes les règles précédentes peuvent être vues comme un cas particulier 

de celle-ci, la difficulté étant le choix de poids w(X). De plus cette règle ne fait pas apparaître 
la gestion fine des conflits partiels telle quelle est considérée dans la règle de Dubois et Prade, 
même s'il est toujours possible d'en tenir compte selon la définition des poids w(X). 

A notre connaissance les seules règles répartissant les conflits partiels sur les éléments dont la 
combinaison crée le conflit, sont la règle du minC proposée par [3] et celles proposées par Dezert 
et Smarandache [T8l 18"]. Ces règles peuvent également être vues comme des cas particuliers de 
la répartition du conflit total de l'équation (fT4l) avec des poids choisis de manière adéquate. La 
règle la plus aboutie de Dezert et Smarandache est la PCR5 donnée pour 2 experts par : 

m PCR5 (X) = m 12 (X) + 

/ mi (X) 2 m 2 (Y) m 2 (X) 2 mi (Y) \ 

^ \ mi (X) + m 2 (Y) + m 2 (X)+m 1 (Y)J ' (15) 

Ye2 e , 

xnY=0 
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où mi2(.) est la combinaison orthogonale pour deux experts, et le dénominateur est non nul. 
L'extension à M experts peut s'écrire : 



?™PCR5W = m Conj (X) 



Aï 



+ E m *w E 



i=i 



H Y ai{k) nx=9 

k — l 



M-1 \ 



E \\{ m ^U)^{j))-T{x=z, mi (x))) 

Y «iU)= z 



(16) 



(y <Tî(1) ,...,y tTi(M _ 1) )e(2 e )«- 1 ^6{x,y CTi(1) ,...,y CTi{M „ 1) } 

où (jj prend des valeurs de 1 à M selon i : 

f <ri(j) = J si j < i, 
\ <r%(j) = j + 1 si j > i, 

et : 



(17) 



T(B,x) = x si .B est vraie, , . 

T(£>, a;) = 1 si i? est fausse. 

Si nous reprenons l'exemple de Zadeh, cette règle donne : 

m PCR5 {A) = 0.486, m PC R5{B) = 0.486, m PC R 5 (C) = 0.028. 

La masse équivalente sur A et B provient des masses initiales des deux experts, pour lever 
l'indétermination il suffit d'une légère différence entre les deux experts. 

3 Une nouvelle règle de combinaison répartissant le conflit 

La règle précédente (équation (fTïïj) ) n'est pas toujours satisfaisante comme nous le montrons 
dans les considérations ci-dessous. Ainsi nous avons proposé une autre extension à M experts 
(équivalente dans le cas de deux experts) : 



mp C R MO (X) = m C on]{X) 

M 



E 



X\m ai Q){Y ai{j) ) 



m 



i=l M-1 

(^(i),..,^( A «))e(2 e ) Af - 1 



M-1 



(19) 



où a est défini par l'équation (fTTI ). 
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Nous pouvons proposer deux règles encore plus générales données par : 

mpcR MOt ( x ) = m c (X) + 

II m ^(i)( y *i(i)) 



M-l 

M 



J2mi(X)f( mi (X)). 



3=1 (20) 



M-l 

i=l 



(y ll ...,KM- 1 )6(2e ) M- lu ^-i } f( mi (X))+ Y.f^U^U))) 



nf = - 1 1 y fe nx=0 j=1 



avec les mêmes notations que dans l'équation (TOI , et / est une fonction croissante définie sur 
]0, 1] et à valeurs dans H + *. 



rapcR MOg ( X ) = m D {X) + 

M-l 

M n m ^u)( Y ^u)) 

^TmiXMmiiX)). —, 

1=1 (Y u ...,Y M . 1 )e(2^^\{X^ } g ( mi (X) + "W^tf)) 

n%-?Y k nx=Q V j=i 

avec les mêmes notations que dans l'équation (fTÔj) . et g est une fonction croissante définie sur 
]0, 1] et à valeurs dans IR + *. 

Par exemple, f(x) = g(x) = x a , avec a G IR + . Sur ces deux dernières règles, la difficulté du 
choix de la fonction / ou g s'apparente à celle du choix des poids dans l'équation ( TT4l) . 

3.1 Remarques sur cette nouvelle règle de combinaison 

- Le problème principal des règles de redistribution du conflit est peut-être la non-associativité. 
Prenons l'exemple de trois experts et de deux classes donné par : 








A 


B 





Expert 1 





1 








Expert 2 








1 





Expert 3 








1 






En fusionnant l'expert 1 et 2 puis l'expert 3, la règle donne : 

m 12 (A) = 0.5, m 12 (B) = 0.5, 

et 

m {12 ) 3 {A) = 0.25, m (12)3 (£) = 0.75. 
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Maintenant, si nous fusionnons l'expert 2 et 3 puis l'expert 1, la règle donne : 

m 23 (A) = 0, m 23 (B) = 1, 

et 

Tn(i2) 3 (A) = 0.5, 771(12)3(5) = 0.5. 

Ainsi le résultat est bien différent. 

De plus avec l'équation (fTïïjl nous obtenons : 

m (m {A) = 1/2, m ( i23)(S) = l/2, 
et avec la règle PCRmo nous obtenons : 

m (123) (A) = l/3, m (123) 0B) = 2/3, 
qui est un résultat plus intuitif. 

L'associativité peut être importante en fusion dynamique. Cependant, lors de la fusion 
dynamique il est courant d'affaiblir les masses afin de tenir compte plus faiblement des 
avis des experts s'étant exprimé il y a longtemps. La PCRmo le fait naturellement, certes 
sans avoir le choix des poids d'affaiblissement, mais il est aussi recommandable dans la 
mesure du possible de procéder à un affaiblissement avant la combinaison. 
Le conflit n'est pas seulement redistribué sur les singletons. En effet, prenons par exemple 
trois experts : 





AU B 


B U C 


AUC 





Expert 1 


0.7 








0.3 


Expert 2 








0.6 


0.4 


Expert 3 





0.5 





0.5 



Le conflit est ici donné par 0.7x0.6x0.5=0.21, avec la règle PCRmo nous obtenons : 

m ( i33)(A) = 0.21, 

m im (B) = 0.14, 
m (123 )(C) = 0.09, 

m ( i23)04 U B) = 0.14 + 0.21.-^ ~ 0.222, 

lo 

m {123) (B U C) = 0.06 + 0.21.— = 0.13, 

18 

m(i23)GA U C) = 0.09 + 0.21.4 - 0.147, 

lo 

m (123) (e) = 0.06. 
Nous obtenons ainsi la classe A quelque soit le critère de décision : 
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A 


B 


AU B 


C 


AuC 


BUC 


e 


«1(123) 





0.21 


0.14 


0.2217 


0.09 


0.1483 


0.13 


0.06 


bel(i2 3 ) 





0.21 


0.14 


0.5717 


0.09 


0.4483 


0.36 


1 


Pl(123) 





0.64 


0.5517 


0.91 


0.4283 


0.86 


0.79 


1 


betP ( i23) 





0.415 


0.3358 


0.7508 


0.2492 


0.6642 


0.585 


1 



- La PCRmo apporte des décisions qui peuvent être différentes des autres règles de combi- 
naison répartissant le conflit. Reprenons l'exemple choisi précédent : 








A 


B 


AUB 


C 


AUC 


BUC 





ITLl 





0.5421 


0.0924 





0.2953 








0.0702 


m 2 





0.2022 


0.0084 





0.6891 








0.1003 


m PCR MO 





0.5421 


0.0924 





0.2953 








0.0702 


bel PC R MO 





0.5421 


0.0924 


0.6345 


0.2953 


0.8374 


0.3877 


1 


PIpCRmo 





0.6123 


0.1626 


0.7047 


0.3655 


0.9076 


0.4579 


1 


betPpcR MO 





0.5655 


0.1158 


0.6813 


0.3187 


0.8842 


0.4345 


1 



Sur cet exemple la classe retenue par la PCRmo sera donc la classe A, contrairement à 
la règle de combinaison de Dubois et Prade (équation (fTBl qui donne la classe C. En 
dehors des considérations philosophiques sur l'intérêt de reporter les conflits partiels sur 
les ignorances partielles puis sur les singletons lors de l'étape de décision ou l'intérêt 
de reporter ceux-ci directement sur les singletons lors de l'étape de combinaison, ce qui 
importe est la règle donnant les meilleures performances. L'exemple précédent montre 
qu'il existe des situations dans lesquelles la décision est différente sans pour autant qu'il 
soit possible de dire quelle règle est la meilleure. Pour ce faire nous allons étudier le 
comportement de ces règles dans le cadre de deux applications où la réalité est supposée 
connue dans la section [4j 

3.2 La combinaison des fonctions de croyance à valeurs réelles en 
répartissant le conflit 

L'article récent de Ph. Smets [22] propose l'extension des fonctions de croyance sur les 
nombres réels. Cette extension ouvre la porte à de nombreuses applications des fonctions de 
croyance dans le domaine continu telles que l'estimation ou le suivi de cibles. Les fonctions de 
masse à valeurs réelles sont définies sur des intervalles, par exemple de IR. Pour les mêmes raisons 
que dans le domaine discret, les experts donnant une masse sur des intervalles disjoints entrent 
en conflit. La combinaison de Smets [22] pour deux experts s'exprimant sur [ai, bi] et sur [02, 62], 
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dont l'intersection non vide est donnée par [a 1; b\] fl [a 2 , b 2 ] = [max(ai, a 2 ), min(&i, b 2 )] = [a, 6] 
est donnée par : 



l*€L r+OO 

«iConjQM]) = / / mi([x,6])m 3 ([o,î/])dyda: 
+ / / m 1 ([a,y])m 2 ([a:,6])d2/da; 

•/ x=— oo J y=b 

ra r+OO 

+mi([a, b]) / / w 2 ([x, y])dydx 

J x=—oo J y=b 



(22) 



+m 2 ([a, 6]) / / m 1 ([x,y])dydx. 

J x=—oo J y=b 

Ainsi le conflit total, étant la masse attribuée sur l'ensemble vide, est donné par : 



mconj(0) = / / / mi([x,z])m 2 ([z,y})dydx 



z=—oo \J x=—oo J y=z 



+ / / m 1 ([z,y})m 2 ([x,z])dydx ) dz. 

' x=— oo J y=z 



(23) 



Cette règle de combinaison reste non idempotente. Ainsi nous pouvons redéfinir l'auto-conflit 
par l'équation Qj. 

La PCRmo peut donc avoir un intérêt, particulièrement dans les cas de conflit important. 
Elle s'écrit dans le cas de deux experts par : 

«ipcrQMD = m Conj ([a,6]) 

oo r+oo f - oc /"+oo 



mi([a,b]) 2 / / m 2 ([x,y])dydx m 2 (M) 2 / / mi([x,y])dydx 

I J x=b J y=x , J x=b J y=x 

' r+oo r+oo ' r+oo r+oo 

mi{[a,b})+ / / m 2 ([x,y])dydx m 2 {[a,b})+ / / mx^x.yYjdydx 

J x=b J y=x J x=b J y=x 

ra ra ra ra 

mi([a,fe]) 2 / / m 2 ([x,y])dydx m 2 ([a,b]) 2 / / mi([x, y])dydx 

i J x=—oo J y=x , J x=—oo J y=x 

~r Ta, Ta t 



(24) 



mi([o,6])+ / / m 2 ([x,y])dydx m 2 ([a,b])+ / / mi([x,y])dydx 

Jx=—ooJy=x J x=— oo J y=x 

3.3 L'espace de définition des fonctions de masse 

Au lieu de considérer l'ensemble de toutes les disjonctions de l'espace de discernement 2 e , il 
est possible de considérer l'ensemble de toutes les disjonctions et de toutes les conjonctions de 
l'espace de discernement, cet ensemble étant noté D e [6]. Nous autorisons ainsi des intersections 
non vide entre deux éléments de l'espace de discernement. Les masses issues des conflits partiels 
apparaissent donc sur des éléments de D @ et ne sont donc plus réparties sur les éléments de 2 e . 
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Selon les applications, il est possible de ne considérer des intersections non nulles seulement 
pour une partie des éléments de 2 e . Ainsi, les règles de combinaison précédemment présentées 
peuvent s'écrire facilement en considérant des classes d'équivalence au lieu des égalités. 

Du fait de la croissance des fonctions de décision (crédibilité, plausiblité et probabilité 
pignistique) la décision ne peut être prise sur les éléments de D @ correspondant à des intersec- 
tions. Cependant, certaines applications, comme en traitement des images [16] où deux classes 
peuvent se rencontrer sur la zone d'intérêt considérée, peuvent nécessiter la manipulation d'un 
tel espace. Dans le cadre des fonctions de masses continues, lorsqu'il est possible de se contenter 
de la combinaison des intervalles de manière discrète, l'intersection des intervalles prend tout 
son sens et là encore il peut être intéressant de considérer D e . 

4 Applications 

Nous présentons ici deux applications dans lesquelles la répartition du conflit peut avoir un 
intérêt. La première est liée à la caractérisation des fonds sous-marin à partir d'images Sonar, 
la seconde concerne la classification de cibles aériennes à partir de données Radar. 

4.1 Classification des images Sonar 

Les images Sonar sont obtenues à partir des mesures temporelles faites en tramant à l'arrière 
d'un bateau un Sonar qui peut être latéral, frontal, ou multifaisceaux. Chaque signal émis est 
réfléchi sur le fond puis reçu sur l'antenne du Sonar avec un décalage et une intensité variable. 
Pour la reconstruction sous forme d'images un grand nombre de données physiques (géométrie 
du dispositif, coordonnées du bateau, mouvements du Sonar, ...) est pris en compte, mais ces 
données sont entachées des bruits de mesures dues à l'instrumentation. A ceci viennent s'ajouter 
des interférences dues à des trajets multiples du signal (sur le fond ou la surface), à des bruits 
de chatoiement, ou encore à la faune et à la flore. Les images sont donc entachées d'un grand 
nombre d'imperfections telles que l'imprécision et l'incertitude. 

Il est alors difficile pour un expert humain de caractériser précisément et avec certitude les 
fonds marins à partir de telles images. Cependant, en milieu sous-marin les capteurs optiques ne 
permettent pas d'imager rapidement et de grandes surfaces les fonds, ce qui est nécessaire pour 
de nombreuses applications telles que la navigation de robots autonomes. Il est donc nécessaire 
de proposer des approches de classifications automatiques des images Sonar. Nous pouvons 
distinguer 5 types de sédiments : les roches, le sable, la vase, les cailloutis (zone de petites 
roches sur du sable ou vase), et les rides (de sable ou de vase, caractérisées par une texture 
particulière). 

La réalité n'étant pas connue précisément, nous avons constitué une base de données à partir 
de l'avis de différents experts en leur demandant d'indiquer la certitude de l'information fournie 
(en terme de sûr, moyennement sûr et pas sûr). La figure Q] présente le résultat de la segmentation 
obtenue par deux experts, les couleurs indiquant les différents types de sédiments pour une 
certitude. Nous constatons ainsi le manque de précision de la segmentation et l'incertitude des 
experts sur leur classification. 
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FlG. 1 - Avis de deux experts sur la même image. 



4.1.1 L'extraction de caractéristiques discriminantes 

La première étape de la classification est l'extraction de paramètres discriminants. Les 
images Sonar étant texturées, nous allons chercher à caractériser cette texture. La texture 
d'une image ou d'une partie d'une image est définie par une ressemblance globale qui peut être 
vue comme la répétition d'un motif plus ou moins grand, pouvant comporter des différences 
locales très importantes. La texture doit donc être caractérisée au niveau de la globalité de la 
partie considérée et non au niveau du pixel. Ainsi, des imagettes de taille variable (par exemple 
16 x 16 pixels, ou 32 x 32 pixels) sont considérées pour la caractérisation de la texture. 

L'extraction de caractéristiques de texture sur des images peut être réalisée à partir de 
différentes méthodes de décomposition de l'image telle que les matrices de co-occurrence, les 
décompositions en ondelettes, par exemple à partir des filtres de Gabor, ou encore les longueurs 
de plage [H]. Une fois l'image décomposée différents paramètres peuvent être calculés tels que 
énergie, entropie, homogénéité, etc. 

Dans cet article nous considérons une décomposition à partir des matrices de co-occurence 
qui par ailleurs ont donné de bons résultats [H]. Une matrice de co-occurrence est obtenue 
par une estimation des probabilités de transition de niveaux gris entre deux pixels situés à une 
distance donnée (ici de un pixel) selon différentes directions d : 0, 45, 90 et 135 degrés. Cette 
matrice Cd est donc une matrice carrée ayant pour taille le nombre de niveaux de gris considéré 
n g . A partir de cette décomposition Haralick [9] a proposé 14 paramètres, nous en avons retenu 
6 : l'homogénéité, le contraste, l'entropie, la corrélation, la directivité, et l'uniformité. 

V homogénéité est donnée par : 

qui a une valeur élevée pour des images uniformes, ou possédant une texture périodique dans 
le sens de la translation d. 

Le contraste qui caractérise des probabilités de transition élevées pour des pixels ayant une 



13 



grande différence de niveau de gris, est donné par 



n.g-1 n a 

^ïH^ E C ^3l (26) 
9 k =° i,i=i,|t-i|=fc 

h' entropie qui reste faible s'il y a peu de probabilités de transition élevées dans la matrice, 
est estimée par : 

n g n g 

l-J2J2 C ^3)HC d (t,j))- (27) 

i=l j=l 

La corrélation entre les lignes et les colonnes de la matrice est calculée par : 

1 



o~ r a. 



^2^2(i-m x )(j -m y )C d (i,j) 

i=l j=l 



(28) 



où m x , m y , a x , et a y sont respectivement les moyennes et écart-types des lignes et de colonnes 
de la matrice C d . 

La directivité sera élevée si la texture présente une direction privilégiée dans le sens de la 
translation d : 

Tlg 

J2C d (i,i). (29) 

1=1 

L'uniformité caractérise la proportion d'un même niveau de gris dans l'image : 

"q 

J2C 2 d (i,i). (30) 
i=i 

Après une étude de ces paramètres sur les images Sonar, nous avons constaté que la taille 
des imagettes a peu d'influence sur la séparation des histogrammes selon les classes. De même 
la direction de translation d n'influence pas ou peu la séparation des classes. En effet, seules les 
rides peuvent avoir des directions privilégiées, cependant nous ne cherchons pas à séparer les 
rides verticales des rides obliques, mais simplement les rides des autres classes. Ainsi, nous avons 
choisi de considérer les six paramètres décrits ci-dessus en prenant la moyenne arithmétique sur 
les quatre directions de d. 

La figure [2] montre la difficulté qu'il peut découler de la classification à partir de tels pa- 
ramètres pour les images Sonar. En effet, nous constatons que pour l'entropie seule la vase se 
distingue des autres sédiments, alors que pour la corrélation il y a un fort recouvrement. Cette 
confusion du type de sédiment par un paramètre entraîne un auto-conflit de ce paramètre, ce qui 
va générer du conflit lors de la combinaison. Ces deux figures résument bien le comportement 
des quatre autres paramètres. 
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FlG. 2 - Corrélation et entropie. 



4.1.2 Modélisation des fonctions de masse 

Une première approche possible pour les fonctions de masse est celle proposée par Denceux 

f m p (Ci\x^' k) )(x) = a ip exp (jid 2 (x,x^)) 
\ m p (Q\x^' h ^)(x) = 1 - a ip exp (%d 2 (x 1 

où Ci est la classe associée à x^'^ (i.e. S 1 '^ = Ci), et x^ ,h > sont les k vecteurs d'apprentissage les 
plus proches de la valeur x du paramètre p. ai P est un coefficient d'affaiblissement, et gammai 
un coefficient de normalisation qui peuvent être optimisés |24| . 

Les fonctions de masse des paramètres de texture peuvent aussi être définies à partir 
des histogrammes normalisés obtenus pour un expert. Notons d%(x) le nombre d'imagettes 
de la classe Ci dont le paramètre de texture p donne x. La normalisation proposée est : 
d%(x)/ m.ax x Ç£2 i dj(x)). Ainsi la masse pour un paramètre p et pour une classe Cj est don- 
née par : 



m p (C, 



m, 



d p Ax) 



max^y^cff^)) 

i 

(6) = l-^m p (Q). 



(32) 



4.1.3 Résultats 

Notre base d'images Sonar est composée de 42 images fournies par le GESMA (Groupe 
d'Etudes Sous-Marines de l'Atlantique). Ces images ont été obtenues à partir d'un Sonar latéral 
Klein 5400 avec une résolution de 20 à 30 cm en azimut et 3 cm en range dans des profondeurs 
allant de 15 m à 40 m. Ces images ont été segmentées par différents experts. L'estimation des 
masses pour chaque paramètre et selon le type de sédiment est effectuée à partir de l'information 
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donnée par un expert, tandis que les tests sont réalisés à partir des informations fournies par 
un second expert. 

Nous ne considérons ici que des imagettes de taille 32x32 pixels ne comportant qu'un seul 
type de sédiment afin d'obtenir des matrices de confusion classiques. Nous obtenons ainsi une 
base de 30294 imagettes pour l'expert fournissant l'apprentissage et de 30745 imagettes pour 
le second expert. Cinq types de sédiments sont considérés : roche (10%), cailloutis (5%), ride 
(13%), sable (27%) et vase (45%). 

Le tableau [2] présente les résultats de classification issue de la fusion selon les règles de 
combinaison orthogonale, de Dubois et Prade et de la PCRmo, des fonctions de masse définies 
pour chaque paramètre par l'équation (T3TÎ ) ou par l'équation (1321 ). Il faut prendre garde à ces 
pourcentages qui peuvent paraître faibles. L'évaluation de la classification des images Sonar est 
délicate, pour une évaluation plus fine le lecteur peut se reporter à [15] . 



Tab. 2 - Pourcentage de bonne classification par type de sédiment. 







Roche 


Cailloutis 


Ride 


Sable 


Vase 


Global 


équation (I3Î1 


mConj 


27.0388 


18.4974 


24.4464 


11.3517 


20.0870 


18.78 


m DP 


16.1290 


22.8230 


20.0174 


21.8504 


18.4912 


19.72 


m PCR MO 


14.0993 


19.1235 


16.0226 


26.5912 


25.4023 


22.64 


équation (f32f) 


^Conj 


23.3323 





4.7741 


79.8283 


79.7774 


61.55 


m DP 


23.3645 





4.7741 


79.8165 


79.7774 


61.55 


m PCR MO 


28.6819 





8.1273 


81.0280 


77.2168 


61.66 



Nous constatons que l'approche proposée par l'équation (13T1) donne de faibles pourcentages 
en général mais que les résultats pour les classes faiblement représentées sont meilleurs que 
ceux donnés par l'équation ( |32l) (par exemple pour les cailloutis ou les rides). La comparaison 
des trois approches de combinaison montre que la répartition du conflit apporte de meilleurs 
résultats qui sont significativement meilleurs même dans le cas de l'équation (l32l) étant donné 
le grand nombre d'imagettes. Dans une optique d'amélioration des pourcentages il est possible 
de combiner les masses issues des différentes méthodes de combinaison ou bien de combiner ces 
classifieurs tel que nous l'avons proposé dans [14] . 

4.2 Classification de cibles Radar 
5 Conclusions 

Nous avons présenté une nouvelle règle de combinaison en répartissant les conflits partiels 
sur les éléments le produisant. Les comportements et propriétés de cette règle ont été étudiés 
en comparaison avec les règles existantes. Cette règle peut être appliquée dans de nombreuses 
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situations lorsque le conflit est important ou non, en discret ou en continue, ou encore dans le 
cadre de la fusion statique ou dynamique. Nous avons présenté les résultats de cette nouvelle 
combinaison dans le cadre de deux applications : en imagerie Sonar et en classification de cibles 
Radar. Les taux de classification montrent que la répartition du conflit apporte de meilleurs 
résultats que les règles précédemment introduites dans la littérature. Ces résultats ne garan- 
tissent pas que les performances de la nouvelle règle de combinaison seront les meilleures dans 
toutes les applications, mais montrent l'intérêt de celle-ci dans certaines applications. 
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